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PREFACE 


The  work  reported  herein  was  conducted  by  the  Arnold  Engineer- 
ing Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC), 
under  Program  Elements  62807F  and  62302F.  This  work  was  co- 
sponsored by  the  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL) 
and  the  AEDC.  The  results  were  obtained  by  ARO,  Inc.,  (a  subsidiary 
of  Sverdrup  & Parcel  and  Associates,  Inc.),  contract  operator  of  AEDC, 
AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project  Num- 
bers VF416  and  VF449.  The  manuscript  (ARO  Control  No.  ARO-VKF- 
TR-74-77)  was  submitted  for  publication  on  September  4,  1974. 

The  authors  gratefully  acknowledge  the  assistance  of  Dr.  M. 

Kinslow  who  supplied  a description  of  his  condensation  calculation 
for  this  report. 
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1.0  INTRODUCTION 


1.1  BACKGROUND 

The  production  of  high-speed  nozzle  expansions  of  gases  of  accu- 
rately known  properties  is  desired  for  such  applications  as  flow -field 
simulation  studies  and  investigations  of  gaseous  rate  kinetics.  For 
many  of  these  applications  the  gas  employed  is  either  nitrogen  (N2) 
or  air,  and  the  gas  density  of  the  expansion  is  sufficiently  high  that 
neither  translational  nor  rotational  temperature  relaxation  is  of  im- 
portance; furthermore,  the  reservoir  enthalpy  is  such  that  the  vib- 
rational mode  enthalpy  is  either  frozen  or  negligible.  The  only  re- 
maining free -stream  source  of  anisentropic  behavior  is  that  of  con- 
densation which  is  a familiar  problem  to  those  who  have  been  forced 
to  install  expensive  drying  or  heating  apparatus,  or  both,  to  inhibit 
condensation.  Coincident  with  the  desire  to  produce  flow  fields  of 
specified  parameters  is  the  requirement  of  experimental  verification 
of  the  flow -field  parameters,  and  one  desires  that  the  measurement 
be  local,  nonperturbing,  and  of  sufficient  completeness  to  enable 
adequate  characterization  of  the  flow  field.  For  many  expansions  of 
practical  interest  this  includes  the  requirement  of  not  only  determining 
the  existence  of  condensation  but  also,  if  present,  the  extent  to  which 
the  macroscopic  flow  properties  are  affected. 


1.2  OBJECTIVES  AND  APPROACH 

Many  otherwise  valuable  techniques  are  not  appropriate  for  diag- 
nostics of  this  commonly  encountered  type  of  flow  field;  e,  g.  , fluores- 
cence techniques  are  inapplicable  because  of  collisional  quenching 
effects;  both  N2  and  O2  are  infrared  (IR)  inactive,  and  consequently, 
no  IR  technique  is  applicable;  interferometric  techniques  are  incom- 
plete in  that  important  flow  parameters  must  be  inferred  or  modeled 
for  agreement  with  the  results.  However,  parameters  of  interest  may 
be  obtained  for  homonuclear  species  at  sufficiently  high  density  by  the 
laser  Raman  (RS)  and  Rayleigh  (RyS)  scattering  techniques  which  utilize 
the  radiation  scattered  by  the  gas  species  from  an  incident  laser  beam 
source  to  determine  the  individual  monomer  specie  densities  and  tem- 
peratures as  well  as  information  concerning  condensation.  When  using 
these  techniques  of  measurement,  the  goals  of  this  work  include  the 
characterization  of  homogeneous  N2  condensation  for  a particular  nozzle 


5 


AE  DC-TR-74-1 14 


i 


configuration  and  for  a range  of  reservoir  pressures,  and  this  charac- 
terization was  to  include  both  axial  and  radial  profiles  of  the  gas  den- 
sity, temperature,  and  pressure.  A further  goal  was  the  determination 
of  the  reservoir  pressure  dependence  of  the  supersaturation  that  the  N2 
gas  could  maintain  for  a conical  nozzle  expansion  from  a room  tem- 
perature source  and  to  find  the  spatial  description  of  the  mass  fraction 
of  condensate  once  clustering  began.  Finally,  by  using  the  Rayleigh 
scattering  results,  it  was  desired  to  determine  the  mean  cluster  size  • 
throughout  the  cluster  growth  region  and  to  attempt  to  obtain  estimates 
of  the  Rayleigh  scattering  depolarization  ratio  as  a function  of  the  clus- 
ter size. 


2.0  THEORY 


2.1  FLOW-FIELD  CALCULATION 

The  flow  field  studied  in  this  investigation  was  an  N2  conical  nozzle 
expansion  whose  isentropic  behavior  external  to  the  nozzle  was  deter- 
mined using  a method  of  characteristics  solution  (MOCS)  from  Ref.  1, 
which  was  performed  using  an  IBM  370/155  computer.  The  nozzle  exit 
plane  conditions  were  the  MOCS  input  data  and  were  determined  using 
the  method  of  Whitfield  (Ref.  2)  for  calculating  the  internal  nozzle  flow 
field  with  boundary  layer  corrections.  The  flow  fields  investigated  were 
all  underexpanded,  and  it  was  assumed  that  the  subsonic  portion  of  the 
boundary  layer  was  accelerated  such  that  the  Mach  number  (M)  of  the 
flow  at  the  exit  plane  nozzle  wall  was  1.01;  this  value  of  M is  faired  into 
the  calculated  M distribution  internal  to  the  nozzle  for  the  region  where 
M > 2.  Various  methods  of  connecting  the  two  regions  were  investigated 
with  no  significant  differences  resulting  in  the  downstream  MOCS  results. 
This  assumption  is  to  be  contrasted  with  overexpanded  flow  fields  where 
the  free-stream  ambient  pressure  influences  the  boundary  layer  internal 
to  the  nozzle.  The  exit  plane  static-pressure  distribution  and  the  en- 
tirely supersonic  M distribution  were  then  used  as  the  starting  line  input 
for  the  MOCS  for  each  reservoir  condition  studied,  and  these  results 
are  shown  along  with  the  experimental  data. 

Since  the  expansion  is  assumed  isentropic,  with  constant  specific 
heat  ratio  (7),  the  pressure  (P) -temperature  (T)  relationship  for  the 
expansion  is 

p/Po  = {X/To)y/(y-l)  (1) 
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where  the  reservoir  conditions  are  denoted  by  the  subscript  (o).  The 
expansion  proceeds  diagramatically  in  the  P-T  plane  as  shown  in 
Fig.  1 and  as  described  by  Eq.  (1),  and  for  some  ranges  of  the  reser- 
voir parameters,  the  expansion  locus  crosses  the  P-T  vapor  pressure 
curve  producing  a saturated  flow  condition.  Since  the  condensation 
phenomenon  is  intuitively  expected  to  be  initially  a three-body  recom- 
bination process  of  finite  rate  constant,  the  expansion  proceeds  along 
the  isentrope  resulting  in  a supersaturated,  metastable  gas  specie. 

By  assuming  sufficient  collisions  beyond  the  saturation  point  (s)  of 
Fig.  1,  exothermic  condensation  onset  occurs  and  molecular  cluster- 
ing begins  and  continues  in  the  expansion  until  the  gas  density  becomes 
sufficiently  low  that  cluster  or  polymer  accretion  ceases.  Beyond  this 
cessation  point  of  Fig.  1,  the  expansion  continues  as  a two-phase  mix- 
ture of  essentially  a constant  specific  heat  ratio  (y').  The  isentropic 
MOCS  provides  the  correspondence  between  the  thermodynamic  param- 
eters and  the  axial  position  of  the  flow  field,  and,  based  on  the  MOCS 
results,  all  expansions  studied  in  this  work  were  such  that  the  satura- 
tion point  (s)  occurred  downstream  of  the  nozzle  exit.  Consequently, 
the  locations  of  the  saturation  point  and  the  onset  of  condensation  or 
anisentropy  were  observable. 

As  indicated  in  Fig.  1,  one  obviously  wishes  to  determine  the 
spatial  location  of  the  onset  point  (0)  as  well  as  the  axial  and  radial 
variation  of  the  temperature,  gas  density,  and  cluster  growth  and  size. 

Vapor  Pressure 


Figure  1.  Diagram  of  expansion  process. 
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22.  DIAGNOSTICS  DESCRIPTION 


Focusing  a laser  beam  of  wave  number  ( vQ ) into  a small  volume 
gaseous  sample  results  in  both  elastic  and  inelastic  photon- molecule 
scattering  RyS  and  RS,  respectively  (Ref.  3).  For  the  present  study, 
only  the  Stokes  rotational  Raman  process  is  utilized,  and  the  intensity 
(Ij)  for  this  RS  process  from  molecular  species  of  rotational  quantum 
number  (J),  apart  from  unimportant  geometrical  and  transmission 
constants,  is  given  by 


It  <*  n[i/Q  — 4B0(J  + 3/2)] 4 4(2J  + 1)S j exp  [— J(J  + l)0R/TR]/qR 


(2) 


where  the  scattering  ensemble  of  number  density  (n)  is  assumed  to  be 
characterized  by  a Boltzmann  distribution  of  rotational  temperature 
(Tr)  and  characteristic  rotational  temperature  (0r  = BQhc/k  = 2.  86  K). 
Anharmonic  effects  are  neglected  in  defining  the  rotational  constant 
(B0)  and  h,  c,  and  k are  Planck's  constant,  the  speed  of  light,  and 
Boltzmann's  constant,  respectively.  Sj  is  the  strength  factor  char- 
acterizing the  photon-molecule  scattering  event.  For  the  flow -field 
regions  studied,  Tr/  0r  > 4 and  qg,,  the  molecular  rotational  partition 
function,  may  be  replaced  by  (Tr/o^)  where  the  symmetry  factor 
a = 2 for  homonuclear  diatomic  specie's  and  1 otherwise.  The 
quantity  £ for  homonuclear  species  of  nuclear  spin  I is  [(1+  1)/(2I+  1); 
1/(21  + 1)]  for  (even;odd)  J species  if  (-jL)^JP  is.  +1,  where  A is  the 
nucleon  number  of  the  atomic  constituent  and  Pis  (+1;-1)  for  the  elec- 
tronic eigenfunction  parity  eigenvalue  for  (symmetric;  antisymmetric) 
inversion  operations.  If  (-l)^P  is  -1,  then  even-odd  J ordering  is 
reversed,  and  if  the  molecule  is  heteronudear,  £ = 1 for  all  J. 
Therefore,  Eq.  (2)  can  be  written  as 


enlijiy^o  - 4Bo(J  + 3/2)]4(2J  + l)  £>jl  = £nC  - J (J  + 1)0r/Tr 


(3) 


where  all  unimportant  constants  are  included  in  C.  The  left-hand  side 
of  Eq.  (3),  when  plotted  versus  J(J  + 1),  yields  by  iterative  calculation 
Tr,  Furthermore,  any  systematic  deviation  from  an  assumed 
Boltzmann  distribution  is  easily  observed. 

The  specie  number  density  (n)  is  determined  by  experimentally 
summing  Eq.  (2)  over  all  detected  transitions  and  employing  an  in-situ 
calibration  for  evaluating  the  proportionality  constant.  For  this  work, 
the  maximum  summation  limit  for  the  n determination  was  J = 24. 
Obviously  by  using  measured  values  of  n and  T^,  the  static  pressure 
(P)  is  immediately  obtainable. 
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For  a condensing  flow  field,  the  gas  sample  is  composed  of  mono- 
mers and  molecular  clusters  of  an  unknown  size  distribution.  If  nj,  t*i, 
and  % are  the  number  density,  polarizability,  and  polarizability  asym- 
metry factor,  respectively,  of  the  i-mer  (Ref.  3),  the  parallel  1(11 ) 
and  perpendicular  1(1)  components  of  the  Rayleigh  scattered  intensity 
from  the  condensing  gas  sample  can  be  written  as 

I(|')  = (Kn0/A4)Si=1(r.i'n0)[af  + (4/45)0?]  (4> 


1(1)  = (Kno/A4)2i=1(nVno)(02/l5)  (5) 

where  nQ  is  the  reservoir  number  density,  K is  a proportionality  con- 
stant, and  the  values  of  I(  II  ) and  1(1 ) have  been  normalized  to  the  inci- 
dent laser  beam  intensity.  Figure  2 shows  the  orientation  of  the  polar- 
ized and  depolarized  scattered  intensities. 
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Figure  2.  Experimental  arrangement. 
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If  the  flow  field  were  isentropic  and  uncondensed,  the  scattered  in- 
tensity values  I(  II ) and  1(1 ) are 

I°(||)  = (Kn0/A4Kn,/n0)°  [af  + (4/45)0?] 

I°(l)  = (Kno/A4)(n1/no)°(0?/15) 

yielding  the  density-  and  space -independent  depolarization  ratio 
I°(l  )/I°(  II ) for  the  monomer 


(6) 

(7) 


p(1)  * fij/loaf  (8) 

2 q 

where  the  0^  term  in  the  denominator,  which  is  0(10"°)  according  to 

Ref.  4,  has  been  neglected  relative  to  the  af  term.  Superior  zero  de- 
notes the  isentropic  condition.  Dividing  Eq.  (5)  by  Eq.  (7),  yields 


I(||)  = I(||)/l°<j |)  « (n^nj0)  + 


(nj/nj^S^nVnjKa./aj)2 


(9) 


2 

where  represents  the  summation  for  i > 2 and  the  0^  terms  which 
are  O{10~3)  relative  to  the  leading  terms  are  again  neglected.  By 
assuming  nj/ni°«  1,  Eq.  (9)  is  written  in  terms  of  a scattering  func- 
tion f as 


f = I(|J)  - 1 « SjCnj/njXaj/fltj)2 


(10) 


and  f is  a direct  measure  of  the  existence  of  clusters  within  the  scatter- 
ing volume.  If  n^,  is  the  total  number  density  of  the  flow  and  if  nj  w n^,. 


nj^nj  = Xj 


(ID 


where  the  i-mer  mole  fraction  is  related  to  the  condensate  mole 
fraction  Xc  by 


Xi  = P(0XC  . i > 2 


(12) 


The  term  p(i)  can  be  loosely  interpreted  as  the  probability  for  the 
existence  of  the  i-mer  in  the  scattering  ensemble. 


The  depolarization  ratio  p = 1(1  )/I(  II  ) can  easily  be  written  as  a 
function  of  the  depolarization  ratios  p(i)  of  the  i-mers,  viz. , 

P'P(,,~[  1 + Sj'(n/n|Xp(,)/p<,))(ai,«,)2l/[l  . '„,)(«. '« ^ (13) 
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or 


p 'pO  % [1  - ijfn,  ■n|Kp(,)  p(1^)(a1'fl|)2]/(f  + 1)  (14) 

If  the  condensate  is  a monodisperse  distribution  in  cluster  size  so  that 
p(i)  = 1 for  i = J and  if  one  assumes  additivity  of  polarizability  for  the 
■weakly  bound  clusters  (a^  = io-i),  Eqs.  (10)  and  (14)  become 

f(J)  = J2XC  (15) 

and 

P‘P{1)  = [1  + KJ)(p(J,'>(,))]/tl  + f(J)l  (16) 

respectively.  The  condensate  mass  fraction  g is 

g = JXe/[l  , (J  - 1)XJ  (17) 

A Poisson  distribution  function  (PDF)  of  mean  cluster  size  (J)  is 
assumed  to  characterize  the  clusters,  rather  than  the  monodisperse 
distribution  (Ref.  5). 


p(i)  = J‘  exp  (— J)  'i ! , i > 2 


(18) 


For  J > 20,  p(l)  is  negligible,  and  its  exclusion  in  Eq.  (19)  is 
justifiable.  For  the  PDF  and  a-x  = iaj. 


f « (J2  + J)XC 


(19) 


and  g is  again  given  by  Eq.  (17).  Equation  (14)  still  requires  the 
assumption  of  a monodisperse  distribution  [p(J)  = 1]  for  evaluation. 

The  preceding  derivations  show  that  a combination  of  RS  and  RyS 
for  investigation  of  a condensing  flow  field  together  yield  monomer  or 
gas  number  density,  rotational  temperature,  detection  of  condensation 
onset,  and  if  the  condensate  mass  fraction  is  determinable,  values  of 
the  mean  cluster  size  and  its  mole  fraction  are  measurable.  Finally, 
the  depolarization  ratio  of  the  predominant  cluster  is  obtained,  yield- 
ing information  of  the  spatial  symmetry  changes  of  the  cluster  during 
molecular  accretion. 
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3.0  EXPERIMENTAL  APPARATUS 


3.1  FLOW-FIELD  GENERATION 

A 4-  by  10 -ft  vacuum  chamber  provided  with  liquid  N2  and  20  K 
gaseous  helium  (He)  pumping  was  used  to  provide  a low  pressure 
(10-7  to  10"3  torr)  background  environment  for  the  conical  nozzle  ex- 
pansion studies.  The  nozzle  source  was  mounted  on  a motor-driven, 
three-dimensional  movement  mechanism,  and  the  nozzle  was  of 
1.  04-mm  throat  diameter,  14.  5-deg  half-angle,  and  exit  area  ratio 
of  13.4.  Figure  2 shows  the  experimental  configuration,  and  one  sees 
that  the  flow  is  in  the  x-direction,  the  laser  beam  injection  is  along. 
+y,  and  scattered  radiation  is  observed  in  the  +z  direction.  Theogas 
used  was  prepurified  N2  of  99.  998  percent  purity,  and  two  -250  A par- 
ticulate filters  were  installed  in  the  inlet  line  to  eliminate  particulate 
matter  and  minimize  heterogeneous  condensation.  The  gas  reservoir 
was  instrumented  with  standard,  calibrated  pressure  and  temperature 
gages. 


3.2  LIGHT-SCATTERING  APPARATUS 

The  laser  source  was  a Coherent  Radiation  Model  52 -B  argon  ion 
laser  which  was  operated  at  514.  5 nm  at  1.  5 w and  1.  0 w power  for 
RS  and  RyS  measurements,  respectively.  For  flow  visualization, 
4.0-w  total  laser  output  power  was  used.  The  incident  radiation 
polarization  was  rotated  along  the  x-direction,  expanded,  and  focused 
onto  the  centerline  of  the  chamber.  The  scattered  radiation  was 
collected  by  an  f/2  lens  collimated  and  focused  onto  the  input  of  a 
0.85-m  double  grating  spectrometer.  For  Rayleigh  scattering  mea- 
surements, HN-22  Polaroid  material  was  placed  between  the  lenses, 
and  a polarization  scrambler  was  placed  immediately  in  front  of  the 
spectrometer  entrance  slit.  The  entrance  slit  aperture  setting, 
collection  optics  magnification,  and  beam  focusing  together  resulted 
in  observation  of  a 1.  5 -mm -long  by  30-  to  50-/um  diameter  cylindrical 
scattering  volume. 

The  detector  was  a thermoelectrically  cooled  EMI-9502B  photo- 
multiplier, and  the  output  was  processed  by  an  Ortec  photon  counting 
system  for  either  digital  display  or  strip- chart  recording. 
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Several  laser  beam  input  apertures  were  used  to  reduce  background 
radiation  resulting  from  forward  scatter  off  optical  components  and 
plasma  light.  Laser  and  viewing  dumps  were  provided  for  further  re- 
ducing background  signals,  and  all  optically  accessible  surfaces  were 
either  painted  with  a flat  black  coating  or  covered  with  a black  flocking 
material. 

A more  complete  description  of  the  experimental  apparatus  is 
given  in  Ref.  5. 


4.0  RESULTS 


4.1  AXIAL  PROFILES 

Figures  3 to  5 show  the  experimental  results  of  the  axial  profiles 
of  Tf{/T0,  n/riQ,  and  I'(  || ) for  an  expansion  at  the  reservoir  parameters 
of  T0  = 284  K and  PQ  = 3.  4,  6.  8,  and  10.  2 atm.  I'(  II ) is  the  ratio  of 
the  axial  value  of  I(  ||  ) to  the  value  of  I(  II  ) at  the  reservoir  conditions. 
The  MOCS  predictions  are  also  shown  as  are  the  saturation  points 
xs  = (x/D)s.  The  MOCS  prediction  for  the  axial  profile  of  n/nQ  shows 
a discontinuity  in  slope  at  x of  approximately  25  as  shown  in  Figs.  3 
to  5.  This  is  due  to  expansion  effects  from  the  nozzle  lip.  The  ex- 
perimental n/nQ  data  are  30  to  40  percent  lower  than  predicted  in  the 
vicinity  of  x ~ 25,  indicating  that  the  theoretical  discontinuity  in  slope 
was  not  realized  in  practice.  The  onset  of  condensation  and  subsequent 
cluster  growth  is  manifested  by  a dramatic  increase  of  I'(  II ) relative 
to  the  isentropic  prediction,  and  these  onset  locations  are  shown  as 
xq  = (x/D)g.  It  is  seen  that  Tj^/Tq  results  are  in  good  agreement  with 
the  MOCS  prediction  for  the  PQ  = 3.4  atm  expansion  for  all  x and  also 
for  the  PQ  = 6.  8 and  10.  2 atm  expansions  for  x < xq.  The  increase  in 
Tr/T0,  which  is  sometimes  as  much  as  50  percent,  due  to  the  release 
of  the  heat  of  recombination  is  obvious  in  Figs.  4 and  5.  Also  shown 
in  Fig.  5 is  the  axial  variation  of  the  depolarization  ratio  (p),  and  it 
is  seen  that  the  monomer  value  of  p is  obtained  prior  to  xq  but  rapidly 
decreases  following  onset  during  the  cluster  growth  region,  thereby 
exhibiting  a configuration  transition  during  growth  to  a more  nearly 
spherical  scatterer.  Figure  6 shows  the  axial  variation  of  I' (II  ) for 
the  P0  range  of  from  1.  93  to  10.  2 atm.  Several  features  warrant 
noting:  xq  moves  nearer  xs  as  P0  increases;  the  magnitude  of  con- 
densate RyS  increases  as  PQ  increases;  except  for  the  region  near 
x~  25,  noted  previously,  I'(  II  ) is  in  excellent  agreement  with  the 
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Figure  4.  Axial  variation  of  number  density,  l'(||),  and  temperature 
for  PQ  = 6.8  atm  and  Tc  = 287  K. 
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Figure  6.  Axial  variation  of  l'(||)  for  all  reservoir  pressures 
investigated. 
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4.2  RADIAL  PROFILES 

Radial  profiles  for  n/nQ  and  Tr/T0  are  shown  in  Fig.  7 for 
PQ  = 10.  2 atm  at  x = 12.  5 and  19.  8,  and  it  is  seen  that  for  r/D  > 0, 
n/nQ  is  systematically  lower  than  the  MOCS  prediction  whereas 
the  Tr/T0  results  are  greater  than  predicted  for  r/D  > 4.  Radial 
profiles  of  I'(  II ) are  shown  in  Fig.  8 for  the  same  expansion  for 
three  values  of  ft,  The  data  for  ft  = 17.  35  are  quite  interesting  in 
that  two  scattering  peaks  are  symmetrically  located  off  the  center- 
line.  From  Figs.  5 and  6,  it  is  seen  that  ft  = 17,35  is  less  than 
Xq  for  r/D  = 0.  However,  Fig.  8 shows  that  onset  has  already 
occurred  for  r/D  > 0,  and  the  p data  shown  in  Fig.  8 support  this 
conclusion.  The  elevated  values  of  Tr/T0  for  r/D  > 0 of  Fig.  7 may, 
in  part,  result  from  this  behavior. 


Figure  7.  Radial  variation  of  number  density  and  temperature  at  two  axial 
positions  for  PG  = 10.2  atm,  T„  = 288  K. 
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Figure  8.  Radial  variation  of  l'(|j)  at  three  axial  positions  for 
P0  = 10.2  atm,  and  radial  variation  of  depolorization 
ratio  and  cluster  size  for  x = 17.35  and  P„  = 10.2  atm. 
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Figures  9 to  11  show  photographic  observation  for  condensation  onset 
for  three  values  of  P0,  and  the  nozzle,  the  dark  isentropic  expansion 
zone  and  the  bright  onset  zone  are  clearly  evident  as  is  filamentary 
structure  within  the  condensation  growth  region. 


Figure  10.  Flow  visualization  photograph,  axial  scan,  P0  = 10.2  atm. 
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Figure  11.  Flow  visualization  photograph,  axial  scan,  PD  = 13.6  atm. 


5.0  ANALYSIS  AND  DISCUSSION 


By  using  the  RS  results  for  n and  Tr,  the  axial  variation  of  P/P0» 
as  shown  in  Fig,  12,  was  determined.  A systematic  increase  in  P/PQ 
relative  to  the  MOCS  prediction  is  observed  for  x > 30  for  PG  = 6.  8 
and  10,2  atm  expansions  and  is  attributable  to  condensation  onset.  It 
is  to  be  noted  that  no  such  increase  is  seen  for  the  3.  4-atm  expansion 
because,  as  Fig.  3 shows,  onset  occurs  beyond  the  largest  x point  for 
which  n and  Tr  data  were  acquired.  The  axial  P-Tr  results  are  shown 
in  Fig.  13  in  the  P-T  plane,  and  the  y = 1.4  isentropic  expansion  is 
indicated.  Again,  the  systematic  deviations  for  the  PQ  = 6.  8 and  10,2 
atm  expansions  are  noted.  Obviously,  radial  pressure  profiles  are 
obtained  using  the  radial  variations  of  n and  T^,  but  no  additional  in- 
sight results. 

The  axial  variation  of  the  RyS  function  (f)  as  determined  from 
Figs.  3 to  5 is  shown  in  Fig.  14  as  a function  of  PQ.  The  extrapolation 
of  f to  zero  yields  an  indication  of  the  spatial  location  (xp)  of  conden- 
sate growth  onset.  The  axial  variation  of  f is  found  to  vary  as 


x = exp  (-‘-bf) 
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By  using  the  N2  vapor  pressure  data  of  Hilsenrath,  et  al.  (Ref.  6) 
and  the  MOCS  for  each  PD  investigated,  the  ratio  of  xq/x s was  deter- 
mined as  a function  of  PQ,  and  the  results  are  shown  in  Fig.  15.  The 
approach  of  the  location  of  condensation  onset  toward  the  saturation 
point  as  PQ  increases  is  obvious.  By  defining  (s)q,  the  isentropic 
supers aturation  pressure  ratio,  as  the  ratio  of  the  MOCS  values  of  P 
at  saturation  and  onset  and  by  defining  (s')$,  the  isentropic  degrees  of 
supercooling,  as  the  difference  of  the  MOCS  values  of  T at  saturation 
and  onset.  Fig.  15  shows  that  at  the  lowest  PQ  investigated  (s)$  exceeds 
200  and  approximately  40  K supercooling  is  possible.  These  param- 
eters decrease  to  4 and  35  K,  respectively,  at  PQ  = 10.  2 atm. 


Figure  15.  Variation  of  xfl/x,,  degrees  of  supercooling,  and  supersaturation 
ratio  with  reservoir  pressure. 
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To  complete  the  analysis  of  the  P0  = 10.  2 atm  expansion,  our 
colleague.  Dr.  M.  Kinslow,  has  graciously  provided  us  -with  a theo- 
retical centerline  calculation  of  the  condensing  expansion*.  This 
calculation  assumes  an  inviscid,  one -dimensional  flow  field  with  no 
mass  transfer  across  the  stream  tube  and  no  condensate  velocity  slip. 

The  constituents'  equation  of  state  is  assumed  to  be  the  perfect  gas 
relation,  and  the  condensate  growth  rate  is  described  by  the  gas- 
condensate  interaction;  i.  e. , no  condensate-condensate  interaction  is 
considered.  Further,  a unity  sticking  or  condensation  coefficient  was 
assumed  for  gas -condensate  collisions.  The  calculations  used  the  un- 
saturated nozzle  exit  plane  conditions  as  the  starting  data,  and  the 
computations  were  performed  using  an  SDS  9300  computer.  The  initial 
size  and  number  density  of  condensation  nuclei  were  unknown  and  were 
adjusted  to  attempt  to  reproduce  as  accurately  as  possible  the  experi- 
mental results.  For  the  PQ  = 10.  2 atm  expansion,  by  using  a mole 
fraction  of  10"4  for  the  condensation  nuclei  and  approximately  5 mole- 
cules per  condensate  nucleus,  satisfactory  agreement  between  the 
calculations  and  measurements  results.  Figure  16  shows  the  calculated 
axial  Tr/T0  dependence  assuming  that  the  gas  rotational  and  trans- 
lational temperatures  are  equal.  Also  shown  is  the  predicted  axial 
variation  of  condensate  mass  fraction  (g). 

By  using  the  results  for  g and  f (Fig.  14)  for  the  PQ  = 10.  2 atm  ex- 
pansion, Eqs.  (17)  and  (19)  were  solved  to  yield  the  mean  cluster  size  J 
at  the  axial  positions  of  x = 17,  20,  24,  and  40,  and  these  are  shown  in 
Fig.  17.  By  using  Eq.  (16)  and  the  depolarization  ratio  results  of  Fig.  5, 
the  corresponding  values  of  p'^  were  determined,  and  Fig.  17  shows  the 
results.  Quite  obviously,  upon  condensation  onset  and  growth,  the  mean 
cluster  size  becomes  larger,  and  the  shape  becomes  more  spherical  as 
expected.  From  these  results,  it  is  seen  that  for  1 < J <30  the  cluster 
depolarization  ratio  (pW))  is  required  to  be  a double-valued  function  of  J, 
exhibiting  a maximum  in  this  region.  By  referring  to  Fig.  8 for  the 
radial  profile  of  the  measured  values  of  p at  x = 17  for  PQ  = 10.  2 atm, 
the  values  of  f were  determined.  This  yielded  as  a function  of  r, 
and,  by  assuming  a unique  relationship  between  p(J)  and  J,  a mean 
cluster  size  can  be  ascribed  to  the  radial  coordinate  (r).  These  results 
are  shown  in  Fig.  8,  and  support  is  given  to  the  assertion  that  off-axis 
condensation  onset  precedes  that  along  the  centerline,  since  the  cluster 
size  increases  as  one  proceeds  off-axis,  and  it  is  only  reasonable  to 
assume  that  smaller  clusters  must  exist  prior  to  the  formation  of  the 
larger  clusters. 

*A  detailed  description  of  this  calculation  is  included  in  Appendix  A. 
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Figure  16.  Axial  variation  of  calculated  condensate  mass  fraction 
and  temperature. 
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Figure  17.  Variation  of  mean  cluster  size  J and  specific  cluster 
depolarization  ratio  p(J)  with  axial  distance,  x. 
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This  phenomenon  is  understood  when  one  considers  the  process  of 
expansion  of  the  nozzle  flow  into  a vacuum.  As  previously  mentioned, 
the  existence  of  the  background  vacuum  is  not  known  by  the  centerline 
molecules  until  the  plus  and  minus  expansion  characteristics  intersect 
the  axis  at  x « 25.  However,  the  off-axis  molecules  obviously  see  the 
necessity  for  increased  expansion  at  smaller  projected  axial  positions 
and,  consequently,  saturate  sooner.  The  earlier  saturation  then  leads 
to  earlier  onset  with  respect  to  the  projected  axial  location,  and  one 
finds  larger  clusters  off-axis  at  a given  axial  position.  This  observa- 
tion clearly  emphasizes  that  the  detailed  characteristics  of  the  nozzle 
are  of  importance  for  a complete  flow-field  characterization  of  conden- 
sation processes.  Only  if  saturation,  onset,  and  cessation  of  growth 
all  occur  prior  to  the  non-simple  characteristic  interaction  region 
will  the  results  be  independent  of  the  background  pressure  into  which 
the  nozzle  flow  expands. 


6.0  CONCLUSIONS 


The  results  of  this  study  have  demonstrated  the  utility  of  Raman 
scattering  density  and  temperature  measurements  for  characterizing  a 
nozzle  expansion  flow  field.  Specifically,  the  axial  and  radial  profiles 
of  n and  Tr  (and  P)  chronicled  the  transition  of  an  N2  expansion  from 
an  isentropic  flow  field  to  a condensing,  anisentropic  expansion.  Dis- 
tinct differences  in  n were  observed  relative  to  the  MOCS  for  the  isen- 
tropic or  metastable,  supersaturated  expansion  for  the  centerline  near 
x = 25.  Also  seen  were  deviations  from  the  MOCS  prediction  for  both 
n and  T in  the  plume  boundary  region.  Although  such  differences  may 
reflect  our  capability  in  handling  the  MOCS,  it  is  suggested  as  possible 
that  such  discrepancies  are  common  in  the  current  practice  of  obtain- 
ing MOCS.  The  axial  values  of  n were  seen  to  be  insensitive  to  con- 
densation growth  as  expected,  and  the  increase  in  Tr  upon  existence  of 
significant  mass  fraction  of  condensate  g was  both  obvious  and  indica- 
tive of  the  magnitude  of  g,  which  was  most  useful  for  the  a priori,  pres- 
ently incalculable  anisentropic  process. 

A combination  of  RS  and  RyS  was  shown  to  yield  values  of  mean 
cluster  size  assuming  only  the  additivity  of  polarizability.  It  is  noted 
that  this  particular  assumption  may  well  become  untenable  for  large 
clusters,  but  the  quantitative  definition  of  a large  cluster  is  at  present 
uncertain.  The  RyS  depolarization  measurements  have  produced  esti- 
mates of  the  depolarization  ratios  of  individual  cluster  sizes  and. 
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furthermore,  have  shown  that  off-axis  condensation  has  significantly 
preceded  the  centerline  phenomenon.  Finally,  photographic  observa- 
tions of  the  condensation  onset,  growth,  and  expansion  regions  show  ' 
definitively  the  rapid  occurrence  of  onset  and  growth,  allowing  one  to 
quickly  characterize  onset  locations  as  a function  of  expansion  source 
parameters.  These  observations  also  indicate  large  cluster,  filamen- 
tary stream  lines. 
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APPENDIX  A 

CONDENSATION  THEORY  AND  CALCULATIONS 


A necessary  condition  for  the  condensation  of  a gas  in  an  expansion 
flow  field  is  the  crossing  of  the  P-T  isentrope  and  the  P-T  vapor  pres- 
sure curve  when  plotted  in  the  P-T  plane.  This  is  illustrated  qualita- 
tively in  Fig.  1.  The  straight  line  illustrates  the  log  P-log  T relation- 
ship for  the  reversible  adiabatic  expansion  of  a perfect  gas 

_p  /xy/(>'-1) 

p0~  W 

In  an  expanding  adiabatic  isentropic  flow,  with  reservoir  conditions 
<Po.  T0),  the  P-T  conditions  at  any  point  in  the  flow  field  will  proceed 
along  the  straight  line  from  the  initial  conditions.  At  some  point  in  the 
flow,  designated  as  point  "s"  in  Fig.  1,  the  P-T  conditions  will  cross 
the  P-T  vapor  pressure  curve  which  describes  the  necessary  conditions 
for  the  gaseous  state  to  be  in  equilibrium  with  the  liquid  state.  As  the 
flow  proceeds  from  this  point,  it  is  in  the  supersaturated  state  with  re- 
spect to  its  liquid  phase.  The  state  of  the  gas  is  thus  metastable  with 
respect  to  its  phase  from  this  point  and  may  continue  to  proceed  along 
the  isentrope  uncondensed.  Or,  conversely,  condensation  may  occur  at 
some  degree  of  supersaturation,  designated  as  point  "0’1  in  Fig.  1,  and 
the  condensation  continues  until  the  condensed  and  gaseous  phases  re- 
establish equilibrium.  The  condensation  process  releases  heat,  and 
thus  from  this  point  of  condensation,  the  initial  restrictions  of  the  ex- 
pansion are  violated,  and,  in  the  P-T  diagram,  the  expansion  proceeds 
along  some  other  path  as  shown  in  Fig.  1. 

In  the  past,  most  research  in  condensation  has  been  concerned  with 
the  condensation  of  a low  vapor  pressure  vapor  in  a carrier  gas,  e.  g. , 
water  vapor  in  the  air.  The  carrier  gas  acts  as  a reservoir  to  maintain 
equality  of  temperature  between  the  static  gas,  the  vapor,  and  the  con- 
densed phase.  Most  of  the  theoretical  work  in  the  past  has,  therefore, 
been  based  on  the  isothermal  assumption  (see  the  excellent  chapter  on 
condensation  by  Wegener,  Ref.  7).  For  the  condition  of  a pure  gas, 
most  investigators  have  simply  applied  the  isothermal  results  where 
the  assumptions  are  usually  not  justified. 

In  this  appendix,  the  equations  describing  the  condensation  of  a pure 
gas  are  developed.  This  analysis  is  based  on  the  following  assumptions: 
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1.  The  condensing  flow  within  a given  stream-tube  is  inviscid 
one -dimensional  flow. 

2.  The  flow  within  the  defining  stream  tube  is  adiabatic  with 
no  mass  transfer  across  the  boundary. 

3.  The  gaseous  components  obey  the  perfect  gas  relation. 

4.  The  condensed  phase  exists  in  the  form  of  spherical  drops 
or  particles  whose  properties  are  given  by  the  bulk  liquid 
or  solid  properties. 

5.  The  mean  directed  velocities  of  the  gas  and  the  condensate 
are  equal  in  the  nozzle  expansion  and  plume. 

6.  The  growth  rate  of  condensate  is  determined  by  the  inter- 
action between  the  condensate  and  the  gaseous  phase,  i.  e. , 
the  interaction  between  the  condensate  particles  is  neg- 
lected. 

7.  The  condensate  particles  are  considered  to  be  in  the  free 
molecule  flow  regime  relative  to  the  uncondensed  gas 
phase. 

8.  The  probability  of  a molecule  being  condensed  upon  striking 
a condensate  particle  is  unity.  However,  re -evaporation 
will  be  permitted. 

9.  All  condensate  particles  are  of  equal  size. 

The  total  mass  flux  at  any  location  along  the  stream  tube  is  given  by 

“ = AVoo<Pi  + P2>  (A-l) 

From  the  second  assumption,  mis  a constant,  i.  e. , 9m/9x  = 0.  The 
net  mass  exchange  between  the  gas  phase  and  the  condensate  particles 
can  be  determined  by  the  application  of  assumptions  four  through  nine. 

The  flux  of  mass  to  a condensate  particle  is  given  by  the  free  mole- 
cule mass  flux  equation  and  is 


V2*n  tb 

Considering  a particle  in  equilibrium  with  its  gas  phase  for  which 
the  ambient  temperature  and  pressure  are  Tc  and  PV(TC,  a),  respectively, 
the  mass  flux  to  the  particle  is  given  by 
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' ■ 4i7 P ^Ci'^Ala"  * - ' 

,v™rrc 

Under  equilibrium  conditions,  this  expression  also  gives  the  evapo- 
ration rate  which  is  a function  only  of  condensate  properties.  There- 
fore, this  expression  for  the  evaporation  is  valid  for  arbitrary- 
ambient  conditions. 


The  rate  of  change  of  mass  in  a. particle  is  given,  by  . ■' 

, 2 da 
4770  a T~ 

....  i .. : . . ,c  s dt  , ■ ■ 

■ . ■ ' i 

Equating  this  expression  to  the  difference  between  the  condensation 
and  evaporation  rates  and  simplifying  yield,  .... 


Pv(Tc.a) 


vRT, 


r V-RTJ 


m — , . da 

- = -pc\x  Tx 


(A-2) 


in  which  the  time  derivative  has  been  changed  to  a distance  deriva-  . 
tive  by  introducing,  the  flow  Velocity.  . ‘ . . ? , ... 

The  conservation  of  energy  for  the  total  flow  may  be  expressed  as 

+ ~y)  + P2^v=o  [cp2'I  M + ~Y  - C(Tc*a)  + cPl  r]  (A_3) 

From  the  first  assumption,  the  left-hand  side  of  Eq.  (A-3)  is  indepen- 
dent of  x,  and  the  mean  latent  heat  of  condensation  is  given  by 


L(Tc,a)  = / L(Tr..a)a2da //  a2da 
a " /a 


The  energy  per  unit  time  liberated  by  the  gas  condensing  on  a 
single  condensate  particle  is 

■luP  .«2 

■ - [L(T  .a)  - C_  (T  - T )] 

v/2r:  RT^  c .'  P1  c • °°  • ■ >'  - 

and  that  required  for  the  evaporating  gas  is 

'in  Pv(Tc,'a)a2  L(Tc.a) 
yj~2nlilc 


33 


AE  DC-TR-74-1 14 


The  rate  of  change  of  energy  of  the  particles  may  be  written  as 


C 


Pc 


dT 


c 


dl 


Equating  this  to  the  net  flux  of  energy  to  a particle  and  simplifying 
give 


j=i[Ln>)  - CP1(T„-  TJ] 


n/«t7 


UTC,.)  .&..C 


(A-4) 


where,  as  before,  the  time  derivative  has  been  changed  to  a derivative 
with  respect  to  flow  distance. 

The  conservation  of  momentum  yields 

"A  "dT  = ^ l^i  + p2)AVi]  (A-5) 

There  is  no  net  momentum  exchange  between  the  gas  and  condensate, 
since  it  was  assumed  that  there  is  no  velocity  differential. 

The  third  assumption  permits  writing  the  perfect  gas  relation  for 
the  uncondensed  gas  and  the  "gas11  of  condensed  particles,  respectively, 
as 


Pi  = Pi  n r 


and 


p2  = P2rvz 

The  static  pressure  due  to  all  particles  is  simply  the  sum  of  these  two 
terms  and  is 


Po=  = RT=>i  + P 2,z)  (A-6) 

Based  on  the  fourth  assumption,  a relation  between  the  number  of 
molecules  per  condensate  particle  (z)  and  particle  properties  can  be 
written  as  follows: 


4lr  Pc 


a 


3 


7. 


3m 


(A-7) 
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Equations  (A-l)  to  (A- 7)  are  basic  equations  governing  the  eight 
unknowns  (P,,,,  T,,,,  V^,  p±,  p£,  a,  z,  and  Tc).  It  is,  therefore,  neces- 
sary that  an  additional  equation  be  obtained. 

As  described  earlier,  as  a gas  expands  and  crosses  the  vapor  pres- 
sure curve,  condensation  commences  on  existing  particles  in  the  flow. 
Condensation  can  occur  on  nuclei  of  foreign  particles  (heterogeneous 
nucleation)  or  on  naturally  occurring  clusters  formed  from  the  parent 
gas  (homogeneous  nucleation).  As  an  equilibrium  gas  or  vapor  is  slowly 
cooled,  dimers,  trimers,  and  in  general  N-mers  are  present  in  equilib- 
rium concentrations  which  increase  as  the  temperature  is  lowered.  For 
a finite  expansion  rate,  clusters  are  formed  in  the  gas,  and,  while  their 
concentration  is  less  than  the  equilibrium  value,  the  equilibrium  con- 
centration is  locally  being  approached  at  a rate  determined  by  flow 
properties. 


It  is  beyond  the  scope  of  the  present  theory  to  analytically  deter- 
mine the  concentration  of  these  clusters  which  subsequently  serve  as 
nucleation  centers.  In  lieu  of  this,  it  is  assumed  that  nucleation  centers 
exist  in  the  flow  in  a given  proportion  and  size.  This  proportion  and  size 
will  be  determined  to  best  fit  the  experimental  condensation  results. 


Let  the  number  fraction  of  condensation  nuclei  in  the  uncondensed 
flow  be  denoted  by  e,  in  other  words 


Nj  + n2 


forP  < P V(T  ,a  -»  ~) 


A consequence  of  assumption  six  is  that  the  number  flux  of  condensate 
particles  is  constant  at  any  point  along  the  stream  tube.  In  other  words, 
a condensation  nucleus  and  later  a condensate  particle  cannot  vanish  or 
new  ones  appear,  but  can  only  grow  (or  decay)  at  the  expense  of  the  gas 
phase. 

This  can  be  described  analytically  as 

i^AV^  = constant 

The  relation  between  number  density  and  mass  density  are 
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and 


Using  these  results  produces  the  final  required  equation 

e 7.  m 

P 2 = (A-8) 

Equations  (A-l)  to  (A-8)  are  the  basic  equations  that  describe  the 
condensing  flow.  In  order  to  solve  these  equations,  it  is  necessary  to 
know  the  following: 

1.  Flow  conditions:  m,  h . A(x) 

o 

2.  Condensate  properties:  Pv(Tc,a),  p L(T  ,a),  C 

”c 

3.  Gas  properties:  R,  Cp,  y 

4.  Nucleation  parameters:  «,  z; 

Equations  (A-l)  to  (A-8)  were  numerically  solved  on  an  SDS  9300 
digital  computer.  A brief  description  of  the  numerical  values  used  in 
the  calculation  follows. 

It  was  assumed  that  the  condensation  effects  are  not  coupled  to  the 
flow-field  calculations.  Therefore,  an  effective  area  for  the  one- 
dimensional condensation  calculation  can  be  obtained  from  the  method- 
of- characteristics  calculation.  For  one-dimensional  flow,  the  absolute 
magnitude  of  the  cross-sectional  area  of  the  stream  tube  is  immaterial, 
and,  for  convenience,  a stream  tube  is  chosen  such  that  the  cross  sec- 
tion at  the  sonic  point  is  unity. 

The  total  mass  flux  in  the  stream  tube  (m)  and  the  total  enthalpy  (hD), 
which  are  independent  of  x,  were  determined  at  a point  where  all  flow 
conditions  were  known.  Upstream  of  the  point  where  the  flow  crosses 
the  vapor  pressure  curve  (Fig.  1),  the  flow  is  isentropic  and  the  condi- 
tions are  identical  to  those  from  the  MOCS  calculations.  The  point 
chosen  to  start  the  condensation  calculation  was  the  nozzle  exit. 

For  the  nitrogen  flows  tested,  the  following  numerical  values  were 
used  in  the  calculations.  For  the  gas  properties. 
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v - 1 * 

R = 2.968  x 106  erg  'gm  K 
Cp,  = 3.KR 

S2  = 3-SR/- 

The  surface  tension,  ct,  in  dyne/ cm  is  given  by 

a = 24.7  - 0.204  T for  T > 63.15  K 

a = H.9  for  T < 63.15  K 

The  latent  heat  of  condensation  in  erg/gm  is  assumed  to  be  given  by 

L = 2.86  X 109  - 1.13  X 107  T - 2'pca  for  T > 63.15  K 

L = 2.79  x 109  - 6.05  x 106  T - 2''pca  for  T < 63.15  K 

The  vapor  pressure  of  condensed  nitrogen  in  torr  is  given  by 
Pv  = exp  (15.859  - 714.67 'T  + 2'pcaRT)  for  T > 63.15  K 

Pv  = exp  (17.635  - 826.83 + 2'pcaRT)  for  T < 63.15  K. 

and  the  bulk  condensate  properties  are  given  by 

C = 1.09  x 10'  erg''gm  K 

Pc 

I 

pc  = 0.873  gm 'cm^  for  T > 63. 15  K 

pt  = 0.947  for  1 < 63.15  K 

While  the  condensation  calculation  could  be  made  along  any  stream 
tube,  the  present  calculations  were  carried  out  only  for  the  centerline 
(r  = 0).  As  stated  earlier,  the  initial  size  and  number  of  condensation 
nuclei  were  considered  as  unknowns.  These  quantities  were  varied  for 
each  flow  condition  in  order  to  obtain  the  best  fit  of  the  calculated  axial 
profiles  of  both  the  Rayleigh  scattering  and  gas  temperature  results  to 
the  experimental  data.  The  values  of  number  fraction  of  nuclei  (e)  and 
initial  number  of  molecules  per  condensate  particle  (z)  obtained  are 
presented  in  Fig.  A-l  as  a function  of  PQ  for  the  nozzle  of  interest. 
Throughout  this  report  these  numerical  results  are  referred  to  as  the 
condensation  calculations. 
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NOMENCLATURE 


A 

B0 


Atomic  number  of  the  atomic  constituent  of  a diatomic 
molecule 

Rotational  constant  for  the  ground  vibrational  level 


b 

C 

c 

D 


Slope  of  the  scattering  function  (f)  when  plotted  versus 
£n(x/  x0) 

A constant  in  Eq.  (3) 

Speed  of  light  in  vacuum 
Nozzle  throat  diameter 


f 


Scattering  function  defined  by  Eq.  (10) 


g 

h 

I 

kid,  ran 

id),  i°(i) 


kid,  ran 


Condensate  mass  fraction 
Planck's  constant 
Nuclear  spin 

Rayleigh  scattered  intensity  from  the  condensing  and 
noncondensing  gas  sample,  respectively,  polarized 
parallel  to  the  incident  beam  polarization 

Rayleigh  scattered  intensity  from  the  condensing  and 
noncondensing  gas  sample,  respectively,  polarized 
perpendicular  to  the  incident  beam  polarization 

Ratio  of  I(  II  ) to  I°(  II  ) and  the  ratio  of  the  axial  value 
of  I(  II  ) to  the  value  of  I(  II  ) at  the  reservoir  conditions, 
respectively 

Intensity  of  the  Jth  Stokes  rotational  Raman  spectral 
line 


J 

K 

k 

M 

MOCS 


Rotational  quantum  number,  or  cluster  size 
A constant  in  Eqs.  (4)  and  (5) 

Boltzmann's  constant 
Mach  number 

Method  of  characteristics  solution 


n,  nQ,  n^.  Number  density  of  N2,  reservoir  number  density  of 

N2,  number  density  of  all  N2  i-mers 
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ni 

■Pj  Po>  Pv 
P 

PDF 

p(i) 

Qr 

RS 

RyS 

A 

r,r 

SJ 

s 

<S>0 

(s)^ 

(s')“ 

T,  Toi  Pr 

Xc 

Xi 

A 

X,  X 
Xg,  X0 


an 


0i 


7,7' 


5 

6r 

X 


Number  density  of  the  i-mer 

Static  pressure,  reservoir  pressure,  vapor  pressure 
of  N2,  respectively 

Electronic  eigenfunction  parity  eigenvalue 
Poisson  distribution  function 
Probability  for  the  existence  of  the  i-mer 
Molecular  rotational  partition  function 
Raman  scattering 
Rayleigh  scattering 

Radial  position  in  the  flow  field,  r = r/D 
Strength  factor 

Saturation  point  in  the  P-T  plane 
Supersaturation  ratio,  <s)g  s Pg/Pv 
Supers aturation  ratio,  (s)q  s Pg/Ps 

o 

Degrees  of  supercooling,  (s'Jg  = Ts  - Tq 

Static  temperature,  reservoir  temperature,  and 
rotational  temperature  of  N2,  respectively 

Condensate  mole  fraction 

i-mer  mole  fraction 

Axial  position  in  the  flow  field,  x = x/D 

Axial  location  of  saturation  and  condensation  onset, 
respectively 

i-mer  polarizability  factor 

i-mer  polarizability  asymmetry  factor 

Ratio  of  specific  heats  before  condensation  and  after 
condensation  ceases,  respectively 

Nuclear  spin  statistical  weight 

Subscript  indicating  condensation  onset  and  the 
characteristic  rotational  temperature,  respectively 

Wavelength  of  scattered  radiation 
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V 


O 


A 

a 

CP 

ho 

L(T,  a) 

m 

• 

m 

N 

P 

PV<TC,  a) 

R 

T 

V 

x 

z 

T 

€ 

P 


Wave  number  of  laser  radiation 
Depolarization  ratio 

Monomer  and  J-mer  depolarization  ratio,  respectively 
Indicates  a summation  for  i = 2,  3,  ... 

Symmetry  factor  for  partition  function 


NOMENCLATURE  FOR  APPENDIX  A 

Cross-sectional  area  of  stream  tube 
Radius  of  condensed  particle 
Specific  heat  at  constant  pressure 
Stagnation  enthalpy 

Latent  heat  of  condensation  at  temperature  (T)  of  a 
drop  of  radius  (a) 

Molecular  mass 

Mass  flux  within  stream  tube 

Number  density 

Gas  static  pressure 

Vapor  pressure  as  a function  of  condensate  tempera- 
ture and  radius 

Gas  constant 

Gas  static  temperature 

Flow  velocity 

Axial  position  in  the  plume  downstream  of  the  nozzle 
throat 

Number  of  molecules  per  condensate  particle 

Ratio  of  specific  heats 

Initial  number  fraction  of  condensation 

Density 
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SUBSCRIPTS 

1,  2 Refer  to  perfect  gas  properties  of  uncondensed  and 

condensed  particles,  respectively 

c Refers  to  condensate  particles 

i Refers  to  initial  condensation  nuclei 

® Refers  to  free -stream  properties 
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